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Abstract 


In  this  paper,  we  investigate  the  applicability  of  a  constraint  satisfaction  problem 
solving  (CSP)  model,  recently  developed  for  deadline  scheduling,  to  more  commonly 
studied  problems  of  schedule  optimization.  Our  hypothesis  is  two-fold:  (1)  that  CSP 
scheduling  techniques  provide  a  basis  for  developing  high-performance  approximate 
solution  procedures  in  optimization  contexts,  and  (2)  that  the  representational  as¬ 
sumptions  underlying  CSP  models  allow  these  procedures  to  naturally  accommodate 
the  idiosyncratic  constraints  that  complicate  most  real-world  applications.  We  focus 
specifically  on  the  objective  criterion  of  makespan  minimization,  which  has  received  the 
most  attention  within  the  job  shop  scheduling  literature.  We  define  an  extended  solu¬ 
tion  procedure  somewhat  unconventionally  by  reformulating  the  makespan  problem  as 
one  of  solving  a  series  of  different  but  related  deadline  scheduling  problems,  and  embed¬ 
ding  a  simple  CSP  procedure  as  the  subproblem  solver.  We  first  present  the  results  of  an 
empirical  evaluation  of  our  procedure  performed  on  a  range  of  previously  studied  bench¬ 
mark  problems.  Our  procedure  is  found  to  provide  strong  cost/performance,  producing 
solutions  competitive  with  those  obtained  using  recently  reported  shifting  bottleneck 
search  procedures  at  reduced  computational  expense.  To  demonstrate  generality,  we 
also  consider  application  of  our  procedure  to  a  more  complicated,  multi-product  hoist 
scheduling  problem.  With  only  minor  adjustments,  our  procedure  is  found  to  signif¬ 
icantly  outperform  previously  published  procedures  for  solving  this  problem  across  a 
range  of  input  assumptions. 
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1  Introduction 


The  problem  of  job  shop  scheduling  to  minimize  makespan  has  been  a  subject  of  exten¬ 
sive  investigation  over  the  years  and  represents  one  of  the  most  well-developed  areas  of 
deterministic  scheduling  theory.  However,  with  few  exceptions,  research  on  job  shop 
makespan  minimization  has  been  restricted  to  solution  of  relatively  small  problems 
under  somewhat  idealized  representational  assumptions.  This  is  due,  in  large  part, 
to  a  methodological  bias  toward  development  of  optimal  solution  procedures.  Our  fo¬ 
cus  in  this  paper  is  on  techniques  for  solving  large  job  shop  problems.  Like  several 
other  recent  efforts  (Adams  et  ah,  1988;  Storer  et  ah,  1992;  Balas  et  ah,  1993),  we 
give  up  the  guarantee  of  optimal  solutions  and  instead  concentrate  on  the  development 
of  efficient  solution  procedures  that  tend  to  optimize  rather  well.  Our  specific  objec¬ 
tive  is  to  explore  the  potential  of  “constraint-posting”  scheduling  frameworks,  which 
have  originated  from  research  in  constraint  satisfaction  problem-solving  (CSP)  in  the 
field  of  Artificial  Intelligence,  as  a  means  for  formulating  and  solving  large  makespan 
minimization  problems.  Generally  speaking,  CSP  models  would  seem  to  offer  an  at¬ 
tractive  approach  to  job  shop  scheduling  problems,  since  they  naturally  accommodate 
more  complex  problem  formulations  (e.g.,  involving  bounded-interval  separation  con¬ 
straints  between  job  steps,  flexible  duration  constraints,  sequence-dependent  setups), 
and  support  a  variety  of  systematic  and  local  search  techniques. 

Job  shop  makespan  minimization  is  a  challenging  problem.  Though  efficient,  poly¬ 
nomial  time  solutions  have  been  found  for  selected,  restrictive  versions  of  the  prob¬ 
lem,  (e.g.,  two  jobs,  two  machines  (Jackson,  1956);  two  machines  and  unit  processing 
times  (Hefetz  and  Adiri,  1982)),  it  is  well  known  that  general  problem  formulations 
are  strongly  NP-hard.  Work  on  optimal  solution  procedures  for  general  makespan 
problems  has  focused  most  heavily  on  the  development  of  implicit  enumeration,  or 
branch  and  bound,  algorithms  (Balas,  1969;  Charlton  and  Death,  1970;  Schrage,  1970; 
Florian  et  ah,  1971;  Lageweg  et  ah,  1977;  Barker  and  McMahon,  1985;  Carlier  and 
Pinson,  1989;  Applegate  and  Cook,  1991;  Brucker  et  ah,  1992).  This  line  of  work  has 
shown  steady  progress.  The  notorious  10-job,  10- machine  problem  originally  posed  by 
(Fisher  and  Thompson,  1963)  was  finally  optimally  solved  by  implicit  enumeration  in 
1989  (Carlier  and  Pinson,  1989),  and  a  more  efficient  branch  and  bound  solution  to 
this  problem  was  subsequently  reported  in  (Brucker  et  ah,  1992).  Yet,  this  improved 
solution  procedure  was  not  able  to  solve  larger  20-job,  10-machine  problems  (running 
for  over  3  days  on  a  Sun  4/20  workstation).  As  observed  by  (Lawler  et  ah,  1989),  the 
applicability  of  implicit  enumeration  schemes  is  limited  to  relatively  small  problems, 
and  their  performance  is  quite  sensitive  to  particular  problem  instances  and  initial 
upper  bound  values. 

Other  research  has  investigated  heuristic  approaches.  Simple  priority  dispatching 
algorithms  are  the  most  representative  and  the  most  widely  used  in  practical  envi¬ 
ronments.  (For  a  general  introduction  to  priority  rules,  see  (Panwalker  and  Iskander, 
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1977).)  While  dispatching  algorithms  are  extremely  fast  and  easy  to  implement,  there 
are  also  drawbacks.  The  performance  of  any  given  rule  is  typically  quite  sensitive  to 
problem  characteristics,  and  it  is  generally  quite  difficult  to  find  one  that  dominates  in 
any  particular  environment.  Such  procedures  are  also  susceptible  to  very  poor  perfor¬ 
mance  in  certain  circumstances,  due  to  the  myopic  nature  of  their  decision-making. 

With  the  rapid  increase  in  computing  power  in  recent  years,  a  growing  body  of 
research  has  focused  on  development  of  more  sophisticated  heuristic  methods,  which 
incorporate  various  forms  of  search  and  aim  at  striking  a  better  cost /performance  trade¬ 
off  than  the  extremes  that  are  provided  by  dispatch  and  optimal  solution  procedures. 
One  notable  approach,  reported  in  (Adams  et  ah,  1988;  Balas  et  ah,  1993),  empha¬ 
sizes  bottleneck  tracking  as  a  heuristic  methodology  for  integrating  optimal  solutions 
to  simpler,  one-machine  subproblems.  A  series  of  shifting  bottleneck  procedures  have 
been  defined  which  have  demonstrated  very  strong  performance  on  a  range  of  previ¬ 
ously  published  benchmarks,  and  provide  a  continuum  of  increasingly  more  accurate 
solution  procedures  at  increasing  computational  expense.  Other  work  has  explored  the 
use  of  various  local  search  techniques  as  a  basis  for  approximate  solution,  including 
simulated  annealing  (Matsuo  et  al.,  1988;  van  Laarhoven  et  ah,  1992),  tabu  search 
(Dell’Amico  and  Trubian,  1991;  Taillard,  1993),  genetic  algorithms  (Della  Croce  et  al., 
1992),  and  general  neighborhood  search  (Storer  et  al.,  1992).  This  work  has  also 
produced  strong  approximate  results  on  previously  studied  benchmarks;  it  has  also  in¬ 
vestigated  and  provided  benchmark  results  for  significantly  larger  makespan  problems 
(Taillard,  1993). 

In  this  paper,  we  propose  a  new  heuristic  procedure  for  solving  job  shop  makespan 
minimization  problems.  Our  approach  is  somewhat  unconventional.  We  start  from  a 
base  procedure  called  TCP  (Precedence  Constraint  Posting)  (Smith  and  Cheng,  1993), 
previously  developed  for  efficient  heuristic  solution  of  job  shop  deadline  scheduling 
problems.  We  reformulate  the  makespan  problem  as  one  of  solving  a  series  of  different 
but  related  common  deadline  problems,  and  define  an  extended  procedure,  Multi-PCP, 
which  employs  PCP  as  a  subproblem  solver.  PCP  has  three  attractive  properties  in  this 
context:  (1)  it  relies  on  simple,  computationally  inexpensive  CSP  heuristics  and  thus 
is  quite  efficient,  (2)  it  performs  only  a  fixed  amount  of  search,  making  its  run  time  a 
predictable  function  of  problem  size  parameters,  and  (3)  the  procedure  can  be  straight¬ 
forwardly  generalized  to  accommodate  more  complex  scheduling  problem  formulations. 
Our  hypothesis  is  that  CSP  scheduling  models  such  as  PCP  can  provide  a  basis  for 
heuristic  makespan  minimization  procedures  that  offer  good  cost/performance  and,  at 
the  same  time,  naturally  extend  to  incorporate  the  more  idiosyncratic  constraints  that 
must  be  enforced  in  many  application  environments. 

To  test  the  viability  of  our  hypothesis,  we  first  analyze  the  performance  of  Multi- 
PCP  on  previously  studied  benchmark  problems,  using  the  above  mentioned  shifting 
bottleneck  procedures  as  a  comparative  base.  On  classical  (small)  benchmark  prob¬ 
lems,  we  find  Multi-PCP  to  perform,  on  average,  better  than  SBl  (the  most  efficient 
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but  least  accurate  shifting  bottleneck  procedure)  and  close  to  SB3  (more  elfective  and 
more  costly)  with  computation  times  equivalent  to  SBl.  Closer  analysis  of  results  in¬ 
dicates  that  the  relative  performance  of  Multi-PCP  depends  significantly  on  the  ratio 
of  number  of  jobs  to  number  of  machines  in  the  input  problem;  for  problems  with  low 
job  to  machine  ratios,  Multi-PCP  is  found  to  consistently  outperform  SBl,  in  several 
cases  outperform  SB3,  and,  on  a  few  problems,  outperform  all  shifting  bottleneck  pro¬ 
cedures.  Conversely,  Multi-PCP  is  found  to  be  less  effective  on  problems  with  high  job 
to  machine  ratios.  These  performance  trends  are  further  confirmed  through  analysis  of 
SBl  and  Multi-PCP  on  more  recently  generated  (larger)  benchmark  problems.  How¬ 
ever,  at  larger  problem  size  levels,  the  sensitivity  of  SBl  to  characteristics  of  specific 
problem  instances  also  becomes  much  more  apparent;  computation  times  of  SBl  are 
found  to  vary  significantly  from  problem  instance  to  instance,  and  the  version  of  SBl 
tested  was  unable  to  solve  all  problem  instances.  The  computational  cost  of  Multi- 
PCP,  in  contrast,  is  seen  to  be  quite  predictable  and  consistently  much  lower  within 
all  problem  size  categories. 

To  demonstrate  the  more  general  applicability  of  Multi-PCP ,  we  also  consider  its  ap¬ 
plication  to  a  more  idiosyncratic  makespan  minimization  problem:  the  multi-product, 
hoist  scheduling  problem  previously  studied  in  (Yih,  1994).  This  problem  is  compli¬ 
cated  by  the  need  to  enforce  bounded  interval  constraints  on  both  operation  processing 
times  and  the  allowable  delay  time  between  consecutive  job  steps,  while  additionally 
accounting  for  the  time  required  to  position  the  hoist  for  any  material  transfer.  We  de¬ 
velop  simple  extensions  to  the  base  PCP  procedure  to  account  for  sequence-dependent 
hoist  travel  times,  and  introduce  a  simple  problem  decomposition  scheme  to  improve 
the  procedure’s  computational  performance.  Over  a  range  of  system  configurations 
defined  by  the  ratio  of  mean  processing  time  to  hoist  speed,  Multi-PCP  is  shown  to 
consistently  and  significantly  outperform  Yih’s  previously  reported  procedure.  More¬ 
over,  this  performance  improvement  is  seen  to  be  invariant  to  the  degree  of  flexibility 
in  the  manufacturing  system,  as  defined  by  the  tightness  of  the  bounded  interval  con¬ 
straints. 

The  remainder  of  the  paper  is  presented  as  follows.  In  the  next  section,  we  first 
consider  formulation  of  job  shop  deadline  scheduling  problem  as  a  CSP  and  summarize 
the  basic  PCP  procedure.  In  Section  3,  we  define  the  extended  Multi-PCP  procedure 
for  makespan  minimization.  The  results  of  our  experimental  study  with  classical  job 
shop  scheduling  benchmarks  are  presented  in  Section  4.  In  Section  5,  we  consider 
application  of  Multi-PCP  to  the  hoist  scheduling  problem.  Finally,  in  Section  6,  we 
summarize  major  points  and  conclusions. 
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1.  Apply  constraint  propagation  to  establish  the  current  set  Vd  of  feasible  values  for 
each  unassigned  variable  d; 

2.  Vd  =  (f)  for  variable  d,  backtrack. 

3.  If  no  unassigned  variables,  or  no  consistent  assignments 
for  all  variables,  quit.  Otherwise, 

4.  Select  an  unassigned  variable  d  to  assign. 

5.  Select  a  value  from  Vd  to  assign  to  d. 

6.  Go  to  step  1  . 


Figure  1:  Basic  CSP  Search  Procedure 

2  CSP  Scheduling  Models  and  the  Job  Shop  Dead¬ 
line  Problem 


Constraint  satisfaction  problem  solving  (CSP)  has  long  been  an  area  of  active  research 
within  the  field  of  Artificial  Intelligence,  and  increasingly,  CSP  models  and  heuristics 
have  been  investigated  as  a  means  for  solving  scheduling  problems  (Cheng  and  Smith, 
1994;  Minton  et  ah,  1992;  Muscettola,  1993;  Sadeh,  1991;  Smith  and  Cheng,  1993). 
Generally  speaking,  a  constraint  satisfaction  problem  (CSP)  is  formulated  as  a  triple 
{V^  D,  C},  where  V  is  a  set  of  decision  variables,  D  a  set  of  domains  for  the  variables 
in  V",  and  C  a  set  of  constraints  on  two  or  more  variables  in  V.  Basic  CSP  solution 
procedures  construct  solutions  through  depth  first  extension  of  partial  assignments, 
according  to  the  following  basic  search  procedure  given  in  Figure  1.  Within  this  search 
procedure,  step  1  is  often  referred  to  as  consistency  enforcement^  and  steps  3  and  4 
are  generally  referred  to  as  variable  ordering  and  value  ordering  respectively.  Specific 
CSP  algorithms  vary  in  the  type  and  level  of  consistency  enforcement  that  is  em¬ 
ployed  (Mackworth,  1977;  Haralick  and  Elliott,  1980),  in  the  mechanisms  incorporated 
for  recovering  from  inconsistent  search  states  (e.g.,  backjumping  (Gaschnig,  1979), 
dependency-directed  backtracking  (Stallman  and  Sussman,  1977),  dynamic  backtrack¬ 
ing  (Ginsberg,  1994)),  and  in  the  heuristics  utilized  for  variable  and  value  ordering. 

CSP  scheduling  research  has  generally  emphasized  development  of  more  specialized 
algorithmic  components,  which  take  advantage  of  the  structure  of  this  particular  class 
of  problems.  Much  of  this  work  has  focused  on  variations  of  the  job  shop  deadline 
problem.  A  job  shop  deadline  problem  involves  synchronization  of  the  production  of  n 
jobs  in  a  facility  with  m  machines,  where  (1)  each  job  j  requires  execution  of  a  sequence 
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of  operations  within  a  time  interval  specified  by  its  ready  time  rj  and  deadline  dj,  and 
(2)  each  operation  Oi  requires  exclusive  use  of  a  designated  machine  Mi  for  a  specified 
amount  of  processing  time  pi.  The  objective  is  to  determine  a  schedule  for  production 
that  satisfies  all  temporal  and  resource  capacity  constraints.  The  job  shop  deadline 
problem  is  known  to  be  NP-Complete  (Garey  and  Johnson,  1979). 

There  are  different  ways  to  formulate  this  problem  as  a  CSP.  Most  frequently,  it 
has  been  formulated  as  a  problem  of  finding  a  consistent  set  of  start  times  for  each 
operation  of  each  job.  Under  this  formulation,  there  are  n  x  m  decision  variables, 
sti,  whose  possible  values  are  sets  of  start  times.  Most  work  with  this  CSP  model 
has  concentrated  on  techniques  for  exploiting  resource  capacity  analysis  as  a  means 
for  restricting  and  directing  search.  ’’Look-ahead”  heuristics  for  variable  and  value 
ordering,  based  on  repeated  computation  of  expected  resource  demand  over  time  and 
identification  of  bottleneck  intervals,  have  been  shown  to  yield  significant  performance 
improvements  over  general  CSP  heuristics  (Sadeh,  1991).  Various  contention  analysis 
techniques  have  also  been  used  to  enhance  consistency  checking  and  early  search  space 
pruning  capabilities  (Sadeh,  1991),  and  to  improve  backtracking  performance  (Xiong 
et  al.,  1992).  Other  work  has  investigated  the  use  of  “iterative  repair”  search  procedures 
(in  contrast  to  the  basic  constructive  search  framework  of  Figure  1),  which  start  with 
an  infeasible  initial  solution  and  attempt  to  incrementally  eliminate  conflicts(Minton 
et  ah,  1992;  Zweben  et  ah,  1990).  However,  such  “repair”  techniques  have  been  shown 
to  perform  rather  poorly  in  comparison  to  constructive  approaches  (Muscettola,  1993). 

The  PCP  scheduling  framework  (Smith  and  Cheng,  1993;  Cheng  and  Smith,  1994) 
and  other  recent  work  in  CSP  scheduling  (Aarts  and  Smith,  1994;  Boddy  and  Goldman, 
1994;  Harvey,  1994;  Muscettola,  1993)  alternatively  start  from  a  problem  representation 
akin  to  a  disjunctive  graph  formulation  (Balas,  1969).  The  problem  is  instead  assumed 
to  be  one  of  establishing  sequencing  constraints  between  those  operations  contending 
for  the  same  resource.  We  define  a  set  of  decision  variables  Orderingij  for  each  {Oi,  Oj) 
such  that  Mi  =  Mj,  which  can  take  on  two  possible  values:  Oi  ■<  Oj  or  Oj  -<  Oi.  In  this 
case,  the  search  proceeds  by  incrementally  “posting”  new  precedence  relations  into  an 
underlying  temporal  constraint  graph  and  propagating  the  consequences  of  each  new 
constraint  to  verify  consistency. 

There  appear  to  be  several  pragmatic  advantages  to  this  second  approach.  By  defer¬ 
ring  commitment  on  specific  start  times  until  they  are  forced  by  problem  constraints, 
a  larger  set  of  possible  extensions  is  retained,  reducing  the  likelihood  of  arriving  at 
an  inconsistent  state.  Likewise,  from  a  solution  robustness  perspective,  precise  start 
time  decisions  are  delayed  (if  possible)  until  needed  at  execution  time.  Finally,  for¬ 
mulation  as  an  ordering  problem  provides  a  more  convenient  search  space  in  which 
to  operate.  The  basic  insight  underlying  PCP  is  that  the  search  benefits  provided  by 
look-ahead  analysis  of  resource  contention  over  time  can  be  obtained  through  much 
simpler,  local  analysis  of  sequencing  flexibility.  In  (Smith  and  Cheng,  1993),  a  configu¬ 
ration  of  variable  and  value  ordering  heuristics  based  on  measures  of  temporal  slack  are 
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shown  to  yield  very  competitive  problem  solving  performance  to  currently  dominant 
contention-based  approaches  (Sadeh,  1991;  Muscettola,  1993)  at  a  small  fraction  of  the 
computational  cost. 


2.1  Problem  Representation 

The  PCP  scheduling  model  can  be  formalized  more  precisely  as  a  type  of  general 
temporal  constraint  network  (GTCN)  (Meiri,  1991).  In  brief,  a  GTCN  T  consists  of 
a  set  of  variables  {Xi,  with  continuous  domains,  and  a  set  of  unary  or  binary 

constraints.  Each  variable  represents  a  specific  temporal  object,  either  a  time  point 
(e.g.,  a  start  time  sti  or  an  end  time  eti)  or  an  interval  (e.g.  an  operation  Oi).  A 
constraint  C  may  be  qualitative  or  metric. 

A  qualitative  constraint  C  is  represented  by  a  disjunction  {Xi  qi  Xj)y  ...\/{Xi  qk  Xj), 
alternatively  expressed  as  a  relation  set  Xi  {qi,  ...,qk}  Xj,  where  qi  represents  a  basic 
qualitative  constraint.  Three  types  of  basic  qualitative  constraints  are  allowed: 

1.  interval  to  interval  constraints  -  The  GTCN  definition  of  (Meiri,  1991)  includes 
Allen’s  13  basic  temporal  relations  (Allen,  1983):  before,  after,  meets,  met-by, 
overlaps,  overlapped-by,  during,  contains,  starts,  started-by,  finishes,  finished-by, 
and  equal.  For  convenience,  we  additionally  include  the  relations  before- or-meets 
and  after- or-met-by,  which  represent  the  union  of  relation  pairs  {before,  meets) 
and  {after,  met-by)  respectively  (Bell,  1989). 

2.  point  to  point  constraints  -  The  relations  identified  in  (Vilain  and  Kautz,  1986), 
denoted  by  the  set  {<,=,>},  are  allowable  here. 

3.  point  to  interval  or  interval  to  point  constraints  -  In  this  case,  the  10  relations 
defined  in  (Ladkin  and  Maddux,  1989)  are  specifiable,  including  before,  starts, 
during,  finishes,  after,  and  their  inverses. 


A  metric  constraint  C  is  represented  by  a  set  of  intervals  {A, ...,  h}  =  {[^i,  6i], ...,  [ak,bk]}. 
Two  types  of  metric  constraints  are  specifiable.  A  unary  constraint  C;  on  point  Xi  re¬ 
stricts  ATi’s  domain  to  a  given  set  of  intervals,  i.e.  {Xi  €  h)  V  ...  V  {Xi  G  h)-  A  binary 
constraint  Cij  between  points  Xi  and  Xj  restricts  the  feasible  values  for  the  distance 
Xj  —  Xi,  i.e.,  {Xj  —  Xi  £  h)  V  ...  V  {Xj  —  Xi  £  h).  A  special  time  point  Xq  can  be 
introduced  to  represent  the  “origin”.  Since  all  times  are  relative  to  JAq,  each  unary 
constraint  Ci  can  be  treated  as  a  binary  constraint  Ca.^ 

^In  Section  5,  we  extend  the  above  definition  of  qualitative,  interval  to  interval  constraints  to 
incorporate  metric  quantifiers,  which  is  necessary  to  model  sequence-dependent  setup  times.  We 
ignore  this  complication  for  now  to  simplify  presentation  of  the  basic  PCP  procedure. 
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A  GTCN  forms  a  directed  constraint  graph,  where  nodes  represent  variables,  and  a 
edge  i  — >•  j  indicates  that  a  constraint  Cij  between  variables  Xi  and  Xj  is  specified. 
We  say  a  tuple  X  =  {xi,...,Xn)  is  a  solution  if  X  satisfies  all  qualitative  and  metric 
constraints.  A  network  is  consistent  if  there  exists  at  least  one  solution.  Figure  2  depicts 
the  constraint  graph  for  a  simple  2  job,  2  machine  deadline  scheduling  problem. 

An  enumerative  scheme  for  solving  a  GTCN  is  given  in  (Meiri,  1991).  Let  a  labeling 
of  a  general  temporal  constraint  network,  T,  be  a  selection  of  a  single  disjunct  (rela¬ 
tion  or  interval)  from  each  constraint  specified  in  T.  In  the  graph  of  Figure  2  there 
are  4  possible  labelings,  owing  to  the  {before- or-meets,  after-or-met-by}  relation  sets 
introduced  to  avoid  resource  contention  between  operation  pairs  (Oi,  02)  and  (O3,  O4). 
Since  any  basic  qualitative  constraint  can  be  translated  into  at  most  four  metric  con¬ 
straints  (Kautz  and  Ladkin,  1991)  (e.g.,  Oi  before- or-meets  Oj  translates  to  eti  <  stj), 
any  labeling  of  T  defines  a  Simple  Temporal  Problem  (STP)  network  -  a  metric  network 
containing  only  single  interval  constraints  (Dechter  et  ah,  1991).  T  will  be  consistent 
if  and  only  if  there  exists  a  labeling  whose  associated  STP  is  consistent. 

For  any  STP  network,  we  can  define  a  directed  edge-weighted  graph  of  time  points, 
Gd,  called  a  distance  graph.  An  STP  is  consistent  if  and  only  if  the  corresponding 
distance-graph  Gd  has  no  negative  weight  cycles.  The  minimal  network  of  the  STP  can 
be  specified  by  a  complete  directed  graph,  called  the  d-graph,  where  each  edge,  ^  -f  j, 
is  labeled  by  the  shortest  path  length,  spij,  from  point  i  to  point  j  in  Gd  (Dechter  et  ah, 
1991).  An  STP  network  can  be  solved  in  0{n^)  time  by  the  Floyd-Warshall’s  all-pairs 
shortest-paths  algorithm,  where  n  is  the  number  of  variables  in  the  STP  network. 
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Thus,  a  simple,  complete  procedure  for  solving  a  GTCN  is  to  enumerate  all  labelings, 
solve  each  corresponding  STP  and  combine  results.  We  can  increase  the  efficiency 
of  this  enumeration  procedure  by  running  a  backtracking  search  over  a  meta-CSP 
network,  whose  variables  correspond  to  arcs  in  the  GTCN  that  can  be  labeled  in  more 
than  one  way  and  whose  domains  are  simply  the  set  of  possible  labelings.  In  the 
case  of  the  deadline  scheduling  problem,  this  leads  to  the  set  of  decision  variables 
V  =  {Orderingij}  previously  identified,  and  a  worst  case  complexity  of 


2.2  The  PCP  Procedure 

The  PCP  scheduling  model  (Smith  and  Cheng,  1993;  Cheng  and  Smith,  1994)  augments 
this  basic  backtracking  search  procedure  to  incorporate  simple  analysis  of  the  temporal 
flexibility  associated  with  each  sequencing  decision  that  must  be  made.  This  analysis 
is  utilized  in  two  ways: 

1.  to  specify  dominance  conditions  that  allow  identification  of  unconditional  deci¬ 
sions  and  early  search  space  pruning,  and 

2.  to  provide  heuristic  guidance  for  variable  and  value  ordering. 


Each  of  these  extensions  is  summarized  below. 

Specification  and  use  of  dominance  conditions  in  PCP  derives  directly  from  the 
concept  of  Constraint-Based  Analysis  (CBA)  originally  developed  in  (Erschler  et  ah, 
1976;  Erschler  et  ah,  1980).  This  work  utilized  calculations  of  the  temporal  slack 
associated  with  an  unordered  operation  pair  to  distinguish  among  cases  where  neither 
ordering  alternative,  just  one  ordering  alternative,  or  either  alternative  remains  feasible. 
For  example,  if  slack{Oi  Oj)  =  I ftj’-esti  —  {pi+pj)  <  0  then  Oi  cannot  be  sequenced 
before  Oj.  These  conditions  are  applied  to  detect  and  post  any  “forced”  sequencing 
constraints  at  each  step  of  the  search,  and  to  detect  inconsistent  solution  states. 

In  (Cheng  and  Smith,  1994),  these  dominance  conditions  are  generalized  to  account 
for  the  wider  range  of  constraints  that  are  specifiable  in  a  GTCN.  Suppose  Orderingij 
is  a  currently  unassigned  variable  in  the  meta-CSP  network,  and  consider  the  d-graph 
associated  with  the  current  partial  solution.  Let  and  Cj  be  the  start  and  end 

points  respectively  of  operations  Oi  and  Oj,  and  further  assume  spij  is  the  shortest 
path  length  from  to  Sj  and  spji  is  the  shortest  path  length  from  tj  to  Si.  Then,  four 
mutually  exclusive  cases  can  be  identified: 


Case  1.  If  spij  >  0  and  spji  <  0,  then  Oi  -<  Oj  must  be  selected. 
Case  2.  If  spji  >  0  and  spij  <  0,  then  Oj  -<  Oi  must  be  selected. 


Case  3.  If  spji  <  0  and  spij  <  0,  then  the  partial  solution  is  inconsistent. 

Case  4.  If  spji  >  0  and  spij  >  0,  then  either  ordering  relation  is  still  possible. 

We  note  that  the  “slack-based”  dominance  conditions  of  (Erschler  et  ah,  1976) 
represent  a  special  case  of  the  above  conditions;  under  classical  job  shop  scheduling 
assumptions  (i.e.,  fixed  processing  times,  simple  job  precedence  constraints)  slack{Oi  -< 

Oj)  =  ^pij-  However,  many  practical  scheduling  problems  require  satisfaction  of  more 
complex  temporal  constraints.  For  instance,  manufacturing  processes  sometimes  place 
limits  on  the  amount  of  time  that  can  elapse  between  consecutive  job  steps  -  e.g., 
if  metal  is  heated  for  subsequent  shaping,  then  shaping  must  occur  before  the  metal 
cools.  Similarly,  a  chemical  bath  operation  may  not  necessitate  a  fixed  amount  of 
time  but  rather  require  a  minimum  processing  time  to  be  productive  and  a  maximum 
time  to  avoid  damage.  Under  such  more  complex  modeling  assumptions,  shortest  path 
information  provides  stronger  dominance  criteria. 

The  second  distinguishing  aspect  of  PCP  is  its  use  of  sequencing  flexibility  analysis 
for  variable  and  value  ordering,  which  dictates  how  the  search  should  proceed  in  the 
undecided  states  (case  4  above).  Intuitively,  in  situations  where  several  Orderingij 
decisions  remain  to  be  made,  each  with  both  possibilities  still  open,  we  would  like  to 
focus  attention  on  the  decision  that  has  the  least  amount  of  sequencing  flexibility.  Since 
the  posting  of  any  precedence  constraint  is  only  likely  to  further  reduce  possibilities 
for  sequencing  other  operation  pairs,  delaying  the  currently  most  constrained  decision 
will  only  increase  chances  of  arriving  at  an  infeasible  solution  state.  With  respect  to 
making  the  selected  ordering  decision,  we  intuitively  prefer  the  ordering  relation  that 
leaves  the  search  with  the  most  degrees  of  freedom. 

One  very  simple  estimate  of  the  sequencing  flexibility  associated  with  a  given  Orderingij 
is  the  minimum  shortest  path  length,  Uij  =  min{spij^  spji)^  which  gives  rise  to  a  variable 
ordering  heuristic  that  selects  the  Orderingij  with  the  minimum  ujij.  This  heuristic 
makes  reasonable  sense;  at  each  step,  the  decision  which  is  closest  to  becoming  forced 
is  taken.  However,  its  exclusive  reliance  on  uJij  values  can  lead  to  problems.  Consider 
two  ordering  decisions  Orderingij  with  associated  shortest  path  lengths  spij  3  and 
spji  =  100,  and  Ordering^!  with  spki  —  4  and  spik  =  4.  In  this  case,  there  are  only 
limited  possibilities  for  feasibly  resolving  Ordering^  and  deferring  this  decision  may 
well  eliminate  them,  while  a  feasible  assignment  to  Orderingij  is  not  really  in  any 
jeopardy. 

To  hedge  against  these  situations,  PCP  instead  bases  variable  ordering  decisions 
on  a  slightly  more  complex  notion  of  biased  shortest  path  length.  Specifically,  hspij 
spij / \/S  and  hspji  =  spji / \/S  are  computed,  where  S  =  min{3p^j ,  spji]  j  xm>x{spij ,  spji} 
estimates  the  degree  of  similarity  between  the  two  values  spij  and  spji.  The  se¬ 
quencing  flexibility  associated  with  a  given  decision  Orderingij  is  redefined  to  be 
uj--  =  min{hspij,hspji),  and  the  decision  selected  during  variable  ordering  is  the  deci¬ 
sion  with  the  minimum  The  value  ordering  heuristic  utilized  in  PCP  simply  selects 
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the  ordering  relation  implied  by  max{bspij,bspji),  i.e.  the  sequencing  constraint  that 
retains  the  most  temporal  flexibility  is  posted. 


2.3  More  Efficient,  Approximate  Procedures 

The  dominance  conditions  and  variable/ value  ordering  heuristics  that  distinguish  the 
basic  PCP  procedure  do  not,  of  course,  change  the  exponential  worst  case  behavior 
of  the  backtracking  search  required  to  guarantee  completeness.  Given  our  pragmatic 
interest  in  solving  large  problems,  we  thus  introduce  two  less-costly,  approximate  so¬ 
lution  procedures  for  later  use.  The  first  variant  is  simply  defined  as  a  backtrack-free 
version  of  the  basic  PCP  procedure.  In  particular,  total  reliance  is  placed  on  the  ability 
of  the  search  to  move  directly  to  a  feasible  solution;  if  Case  3  above  is  ever  encountered 
(i.e.,  no  feasible  ordering  for  a  given  ordering  decision),  the  search  simply  terminates 
in  failure  (and  does  not  produce  a  solution).  The  effectiveness  of  this  partial  solu¬ 
tion  procedure,  which  we  will  refer  to  as  “Simple  PCP”  below,  was  demonstrated  in 
(Smith  and  Cheng,  1993)  on  a  set  of  previously  published  CSP  scheduling  benchmark 
problems. 

We  also  define  a  second  variant,  referred  to  below  as  “Simple  PCP  with  Relaxation” 
which  extends  Simple  PCP  in  the  following  manner.  Whenever  an  ordering  decision  is 
recognized  as  Case  3,  the  unresolvable  decision  is  set  aside,  and  the  search  is  allowed 
to  proceed  with  other,  still  resolvable  ordering  decisions.  Once  all  feasibly  resolvable 
decisions  have  been  made,  the  set  U  of  unresolvable  (Case  3)  decisions  is  then  re¬ 
considered.  For  each  Orderingij  in  U,  deadlines  di  and  dj  are  relaxed  (increased)  by 
\max{spij,spji)\  and  the  corresponding  precedence  relation  (which  is  now  feasible)  is 
posted.  This  second  approximate  procedure  thus  always  produces  a  solution,  albeit  one 
that  may  not  satisfy  all  original  problem  constraints.  Both  approximate  procedures 
can  be  seen  to  have  worst  case  time  complexity  of  On^\V\^  where  |R|  is  the  number  of 
ordering  decisions  that  must  be  made. 


3  MULTI-PCP 


While  constraint  satisfaction  scheduling  procedures  such  as  PCP  have  been  elfectively 
applied  to  complex  deadline  scheduling  problems,  their  applicability  to  more  com¬ 
monly  studied  problems  of  schedule  optimization  is  not  obvious.  In  this  section,  we 
focus  specifically  on  the  problem  of  makespan  minimization  and  propose  one  possible 
approach  to  incorporating  these  techniques. 

Our  approach  is  motivated  by  the  concept  of  problem  duality  exploited  in  the  MUL- 
TIFIT  algorithm  (Coffman  et  ah,  1978)  in  the  context  of  multiprocessor  scheduling. 
Suppose  that  we  are  given  an  instance  of  a  makespan  problem,  denoted  by  11^(7)  where 


I  represents  the  problem  data  associated  with  this  problem  instance.  If  we  know  the 
minimum  makespan  for  to  be  then  we  can  reduce  to  a  special 

deadline  problem  !!£)(/,  d),  where  each  job  is  assigned  a  0  ready  time  and  a  common 
deadline  d,  with  d  =  C'^aa;-  ^  —  ^maxi  assured  that  a  feasible  solution  to 

d)  exists.  More  important,  defines  a  unique  common  deadline  such  that  for 
d  <  n£)(/,  d)  has  no  feasible  solution.  This  dual  relationship  between  problems 

and  n£)(/,d)  implies  that  the  makespan  problem  njvf(/)  can  be  reformulated 
as  a  problem  of  finding  the  smallest  common  deadline,  dmim  for  which  Il]:){I,d)  has  a 
feasible  solution. 

Given  an  algorithm  for  optimally  solving  the  deadline  problem  !!£)(/,  d),  it  is  straight¬ 
forward  to  construct  an  search  procedure  for  determining  dmin  (and  its  associated 
schedule).  We  start  with  known  upper  and  lower  bounds  du  and  di  on  the  common 
deadline  dmin',  at  each  step,  we  attempt  to  solve  !!£)(/, d)  for  d  =  {djj  +  di)/2.  If  a 
feasible  schedule  is  found,  du  becomes  d;  otherwise,  du  becomes  d.  We  continue  the 
search  until  du  =  du,  retaining  the  schedule  with  the  best  makespan  as  we  go. 

There  is  a  complication,  however,  in  utilizing  this  binary  search  procedure  in  con¬ 
junction  with  a  heuristic  deadline  scheduling  procedure.  In  particular,  the  search  may 
fail  to  yield  the  best  solution  if  the  deadline  scheduling  procedure  does  not  ensure 
monotonicity  in  solution  results  across  an  interval  of  common  deadlines.  This  property 
implies  that  if  a  feasible  solution  cannot  be  found  for  a  given  common  deadline  di, 
then  a  solution  will  also  not  be  found  for  any  common  deadline  d2  <  di,  and  like¬ 
wise  if  a  solution  is  found  for  a  given  di,  then  a  solution  will  also  be  found  for  any 
d2  >  di.  It  is  not  difficult  to  construct  examples  which  demonstrate  that  neither  of 
the  simple,  one-pass  PCP  procedures  defined  in  Section  2.3  possess  this  property,  and 
consequently  the  assumptions  underlying  use  of  binary  search  are  no  longer  valid.  For 
this  reason,  we  instead  define  our  extended  makespan  minimization  procedure  in  terms 
of  a  more  conventional  k-iteration  search;  the  approximate  PCP  procedure  (either  vari¬ 
ant)  is  applied  k  times  with  different  common  deadlines  evenly  distributed  between  du 
and  du-  While  k-iteration  search  obviously  also  provides  no  guarantee  of  finding  the 
optimal  solution,  empirical  analysis  has  indicated  that,  with  proper  selection  of  k,  use 
of  k-iteration  search  leads  to  consistently  better  makespan  minimization  performance. 

The  only  remaining  issue  concerns  initial  establishment  of  upper  and  lower  bounds 
on  dmin-  A  lower  bound  du  is  provided  by  the  procedure  originally  described  in  (Flo- 
rian  et  ah,  1971),  where  each  machine  is  sequenced  independently  in  order  of  earliest 
operation  start  times  and  the  maximum  job  completion  time  is  then  selected.  An  up¬ 
per  bound  du  can  be  obtained  through  application  of  one  or  more  priority  dispatch 
rules.  In  the  experiments  reported  below,  a  set  of  six  priority  rules  -  SPT,  LPT,  LFT, 
EFT,  MOR,  and  LOR  -  were  applied,  taking  the  best  makespan  generated  as  du-  The 
complete  algorithm,  referred  to  as  MULTI-PCP,  is  given  in  Figure  3. 
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1.  Compute  upper  bound  du  and  lower  bound  di. 

Set  Best-Makespan  =  du- 

Set  /  =  0. 

2.  If  /  >  k,  stop. 

Otherwise,  set  common  deadline  d  =  du  —  I{du  —  du)lk. 

3.  Apply  PCP  (either  variant)  and  compute  the  makespan  M. 

If  M  <  BestJdakespan,  set  BestJAakespan  =  M. 

If  Best -Makespan  =  du,  stop.  We  have  the  optimal  solution. 
Otherwise,  set  /  =  /  +  1.  and  go  to  2. 


Figure  3:  The  Multi-PCP  procedure 

4  A  Benchmark  Performance  Study 


In  this  section,  we  empirically  analyze  the  performance  of  Multi-PCP  on  two  sets  of 
job  shop  scheduling  benchmark  problems  previously  studied  within  the  literature.  The 
first  problem  set,  referred  to  below  as  the  “small”  problem  set,  consists  of  39  job  shop 
problems  with  sizes  varying  from  6-job  by  6-machine  to  15-job  by  15-machine.  The  first 
three  problems,  Mt06,  MtlO,  and  Mt20,  are  the  long  standing  problems  of  (Fisher  and 
Thompson,  1963).  The  remainder  are  taken  from  the  40  problems  originally  created 
by  (Lawrence,  1984);  of  these  40  problems,  we  include  only  the  36  problems  for  which 
optimal  solutions  have  been  obtained.  The  second  set  of  benchmark  problems,  which 
we  designate  as  the  “large”  problem  set,  are  the  problems  more  recently  defined  by 
(Taillard,  1993).  This  set  consists  of  80  larger  job  shop  problems  with  sizes  ranging 
from  15-job  by  15-machine  to  100-job  by  20-machine.  For  each  problem  in  this  set, 
Taillard  reported  the  “best  solution”  obtained  with  a  tabu  search  procedure  that  was 
run  for  extended  time  intervals. 

We  take  as  a  principal  comparative  base,  the  shifting  bottleneck  family  of  proce¬ 
dures,  SBl,  SB3  and  SB4  (Adams  et  ah,  1988;  Balas  et  ah,  1993),  which  provides  a 
series  of  increasingly  more  accurate  approximate  procedures  for  makespan  minimiza¬ 
tion  at  increasingly  greater  computational  expense.  On  the  small  problem  set,  we 
compare  the  performance  of  each  of  these  procedures  and  Multi-PCP  in  terms  of  two 
measures:  %  deviation  of  generated  solutions  from  the  optimal  solution  and  amount 
of  computation  time  required  for  solution.  Results  for  SBl  were  obtained  on  a  Sun 
SPARC  10  workstation  using  an  implementation  kindly  provided  to  us  by  Applegate 
and  Cook  (for  detail  please  see  (Applegate  and  Cook,  1991)).  Results  for  SB3  and 
SB4  were  taken  from  (Balas  et  ah,  1993),  with  the  reported  Sun  SPARC  330  compu¬ 
tation  times  translated  to  reflect  expected  performance  on  a  SPARC  10.  Multi-PCP 
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Table  1:  %  deviation  from  optimal  solution  for  Multi-PCP,  SBl,  SB3,  and  SB4  across 
small  benchmark  problem  categories _ 


Job  X 

Multi-PCP 

Multi-PCP 
w/  Relax 

SBl 

SB3 

SB4 

Machine 

mean 

a 

mean 

cr 

mean 

a 

mean 

cr 

mean 

(7 

mt06 

0.00 

- 

0.00 

- 

7.27 

- 

0.00 

- 

0.00 

- 

mtlO 

2.04 

- 

2.04 

- 

2.37 

- 

5.48 

- 

1.08 

- 

mt20 

8.76 

- 

2.32 

- 

5.41 

- 

2.92 

- 

2.92 

- 

10  X  5 

1.53 

1.58 

L47 

1.46 

1.59 

1.84 

1.44 

2.05 

1.44 

2.05 

15  X  5 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

20  X  5 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

10  X  10 

2.44 

1.75 

2.44 

1.75 

4.94 

5.32 

3.16 

2.33 

2.25 

1.17 

15  X  10 

4.29 

2.34 

3.78 

1.99 

6.34 

2.60 

2.72 

1.91 

2.72 

1.91 

20  X  10 

3.12 

2.65 

3.12 

2.65 

6.57 

7.26 

1.37 

1.23 

0.96 

1.26 

30  X  10 

0.30 

0.58 

0.30 

0.58 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

15  X  15 

4.15 

0.73 

4.15 

0.73 

6.16 

2.10 

3.09 

0.92 

3.01 

0.97 

All 

1.93 

1.20 

1.72 

1.15 

3.01 

2.39 

1.51 

1.05 

1.24 

0.92 

computation  times  were  also  obtained  on  a  SPARC  10.  All  procedures  considered  were 
implemented  in  C.  For  the  large  problem  set,  we  alternatively  compute  the  %  deviation 
from  the  best  tabu  search  solution  produced  by  (Taillard,  1993)  as  the  performance 
criterion.  Cost/performance  comparison  across  these  problems  is  restricted  to  Multi- 
PCP  and  SBl,  since  results  have  not  been  reported  for  SB3  and  SB4  on  this  problem 
set. 

We  consider  the  performance  of  two  versions  of  the  Multi-PCP  procedure,  defined 
by  incorporating  either  Simple  PCP  or  Simple  PCP  with  Relaxation  as  the  base  CSP 
scheduling  procedure  (see  Section  2.3).  These  two  configurations  are  referred  to  below 
as  “Multi-PCP”  and  “Multi-PCP  with  Relaxation”  respectively.  In  both  configura¬ 
tions,  the  bound  k  on  the  number  of  iterations  performed  was  set  to  8  for  all  runs. 


4.1  Computational  results  on  the  small  benchmark  problems 

Tables  1  and  2  summarize  the  performance  results  obtained  on  the  small  benchmark 
problem  set.  Performance  is  indicated  individually  for  the  3  problems  of  Fisher  and 
Thompson  and  aggregated  according  to  problem  size  for  the  36  problems  of  Lawrence. 
Associated  computation  times  are  given  in  Table  2.  Computation  times  were  found  to 
be  identical  for  both  Multi-PCP  configurations  at  the  level  of  precision  given  in  Table 
2  and  are  thus  listed  only  once.  Detailed  results  for  each  individual  problem  are  given 
in  Tables  5,  6,  and  7  in  Appendix  A. 

Considering  first  the  performance  of  the  two  Multi-PCP  configurations  tested,  we 
see  that  augmenting  the  base  PCP  procedure  to  produce  ’’relaxed”  deadline  solutions 
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10  machines 


Figure  4:  Mean  %  deviation  from  optimal  solution  for  increasing  job/machine  ratio  on 
small  benchmark  problems 

when  feasible  solutions  are  not  found  yielded  improved  solutions  in  only  a  small  num¬ 
ber  of  problems.  In  these  isolated  cases,  however,  the  improvement  provided  by  the 
extended  procedure  was  sometimes  substantial;  for  the  mt20  problem  of  Fisher  and 
Thompson,  %  deviation  from  the  optimum  was  reduced  from  8.76  to  2.32,  matching 
the  best  solution  found  for  this  problem  by  any  of  the  shifting  bottleneck  procedures. 
Since  the  extended  Multi-PCP  with  Relaxation  procedure  incurs  virtually  no  additional 
computational  cost,  we  restrict  attention  to  the  results  obtained  with  this  configuration 
in  our  discussion  below. 

As  shown  in  Table  1,  the  makespan  minimization  performance  of  Multi-PCP  on 
the  small  problem  set  falls  within  the  performance  continuum  defined  by  the  shifting 
bottleneck  procedures.  On  average,  Multi-PCP  is  seen  to  perform  better  than  SBl  and 
very  close  to  SB3,  with  SB4  yielding  the  best  overall  makespan  performance.  Relative 
performance  was  found  to  vary  across  different  problem  subsets.  On  the  three  classic 
Fisher  and  Thompson  problems,  Multi-PCP  found  equivalent  or  better  solutions  than 
both  SBl  and  SB3  in  all  cases,  and  failed  to  match  the  performance  of  SB4  in  just  one 
case.  There  is  little  difference  in  performance  on  the  very  small,  6-machine  problems; 
all  procedures  produce  optimal  or  near  optimal  solutions  in  these  problem  categories. 

Consideration  of  results  on  the  larger,  10-machine  problem  categories  reveals  per¬ 
haps  the  most  significant  comparative  performance  trend.  These  results  are  graphically 
depicted  in  Figure  4  in  terms  of  increasing  ratio  of  number  of  jobs  to  number  of  ma¬ 
chines  in  the  input  problem.  For  problems  with  low  job  to  machine  ratios,  Multi-PCP 
exhibits  its  strongest  comparative  performance.  In  the  case  of  the  10x10  problem 
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Table  2:  CPU  time  (in  seconds)  for  Multi-PCP,  SBl,  SB3,  and  SB4  across  small 
benchmark  problem  categories 


Job  X 

Multi-PCP 

SBl 

SB3 

SB4 

Machine 

mean 

(7 

mean 

(7 

mean 

(T 

mean 

a 

mt06 

0.05 

- 

0.12 

- 

0.78 

- 

1.45 

- 

mtlO 

0.38 

- 

0.72 

- 

2.21 

- 

7.76 

- 

mt20 

1.38 

- 

0.28 

- 

1.99 

- 

3.62 

- 

10  X  5 

0.13 

0.09 

0.12 

0.01 

0.40 

0.20 

0.56 

0.28 

15  X  5 

0.22 

0.26 

0.15 

0.02 

0.12 

0.04 

0.12 

0.04 

20  X  5 

0.07 

0.07 

0.16 

0.03 

0.14 

0.04 

0.14 

0.04 

10  X  10 

0.26 

0.03 

0.67 

0.13 

1.61 

0.16 

3.20 

0.25 

15  X  10 

1.04 

0.10 

1.20 

0.18 

3.41 

0.98 

5.74 

2.52 

20  X  10 

2.51 

0.13 

1.61 

0.26 

3.86 

2.27 

7.50 

5.35 

30  X  10 

4.85 

3.96 

2.58 

0.58 

4.13 

1.76 

4.13 

1.76 

15  X  15 

1.29 

0.17 

4.55 

1.30 

13.34 

0.49 

26.79 

1.05 

Average 

1.19 

0.60 

1.21 

0.31 

2.96 

0.74 

5.29 

1.41 

category,  Multi-PCP  performed  better  on  average  than  both  SBl  and  SB3,  and  very 
close  to  SB4.  Conversely,  Multi-PCP  was  found  to  be  less  effective  (comparatively) 
on  problems  with  high  job  to  machine  ratios.  On  the  30x10  problems  (which  turn  out 
to  be  the  easiest  10-machine  problems  for  all  procedures),  all  three  shifting  bottleneck 
procedures  were  able  to  obtain  optimal  solutions,  whereas  Multi-PCP  failed  to  find  the 
optimum  for  2  of  the  5  problems  in  this  category  (see  Table  6). 

From  a  computational  perspective.  Table  2  shows  that  the  solution  times  achieved 
on  this  problem  set  by  Multi-PCP  are  comparable  overall  to  those  of  SBl.  Table  2 
also  confirms  the  increasingly  higher  computational  demands  of  SB3  and  SB4,  on  the 
15-machine  problems,  the  average  solution  times  for  SB3  and  SB4  are  approximately 
10  and  20  times  larger  (respectively)  than  those  of  Multi-PCP. 


4.2  Computational  results  on  the  large  benchmark  problems 


Table  3  extends  the  performance  comparison  of  Multi-PCP  and  SBl  to  the  larger 
problem  set  of  Taillard.  Again,  performance  results  are  summarized  by  problem  size 
category,  in  this  case  relative  to  the  best  solutions  reported  by  Taillard  s  tabu  search 
procedure.  Figure  5  graphically  depicts  average  solution  performance  by  increasing  job- 
to-machine  ratio  as  before  for  both  15-machine  and  20-machine  problem  categories.  We 
note  that  due  to  excessive  memory  requirements,  the  SBl  implementation  was  able  to 
solve  only  9  of  the  10  problems  in  each  of  the  50x15  and  50x20  problem  sets,  and  was 
unable  to  solve  any  of  the  largest  100x20  problems  on  our  Sun  SPARCIO  workstation. 
The  average  performance  numbers  listed  in  Table  3  and  Figure  5  for  these  categories 
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Table  3:  %  deviation  from  the  best  solution  for  Multi-PCP  and  SBl  across  large 
benchmark  problem  categories _ 


Job  X 
Machine 

Multi-PCP 

Multi-PCP 
w/  Relax 

SBl 

mean 

cr 

mean 

a 

mean 

(T 

15  X  15 

5.74 

0.74 

5.57 

0.78 

9.00 

2.05 

20  X  15 

7.52 

1.82 

7.27 

1.49 

10.15 

2.14 

30  X  15 

9.88 

2.57 

9.65 

2.47 

8.38 

3.06 

50  X  15* 

7.39 

2.31 

7.21 

2.26 

2.66t 

1.57 

20  X  20 

7.60 

1.54 

7.33 

1.33 

9.98 

2.29 

30  X  20 

11.76 

2.43 

11.64 

2.37 

13.05 

2.55 

50  X  20* 

8.77 

0.96 

8.34 

1.06 

5.33t 

1.82 

100  X  20 

4.88 

1.41 

4.88 

1.41 

+ 

- 

All* 

8.38 

1.72 

8.14 

1.65 

8.36 

2.21 

t  SBl  able  to  solve  nine  out  of  ten  problems, 
t  SBl  unable  to  solve  any  of  the  100x20  problems. 

*  Average  performance  is  measured  with  respect  to  problems  solved  by  both  procedures. 


reflect  only  those  problems  that  were  successfully  solved.  Corresponding  average  com¬ 
putation  times  by  problem  category  are  given  in  Table  4.  Detailed  performance  on 
each  individual  problem  is  given  in  Tables  8,  9,  and  10  in  Appendix  A. 

Ignoring  the  scalability  problems  encountered  with  the  tested  SBl  implementation, 
the  results  at  larger  problem  sizes  make  much  more  explicit  the  comparative  perfor¬ 
mance  trends  observed  at  the  10-machine  problem  level.  Multi-PCP  is  seen  to  consis¬ 
tently  outperform  SBl  at  low  job-to-machine  ratios,  while  the  inverse  is  true  at  high 
job-to-machine  ratios.  Both  procedures  achieve  increasingly  better  solutions  at  higher 
job-machine  ratios  (consistent  with  Taillard’s  observation  that  these  problems  are  eas¬ 
ier),  but  in  no  cases  does  either  Multi-PCP  or  SBl  achieve  the  best  solutions  generated 
by  extended  Tabu  search. 

Examination  of  relative  computational  costs  (Table  4)  indicates  some  additional 
scalability  trends  and  tradeoffs.  Multi-PCP  was  found  to  consistently  produce  solutions 
in  less  computation  time  than  SBl;  the  largest  differential  (roughly  4  times  as  fast) 
was  observed  in  the  problem  categories  with  the  smallest  job-to-machine  ratios,  and, 
on  average,  Multi-PCP  obtained  solutions  in  about  half  as  much  CPU  time  as  SBl. 
Multi-PCP  was  also  found  to  be  much  more  predictable  with  respect  to  computational 
cost.  As  shown  in  Table  4,  the  variance  in  Multi-PCP  solution  times  across  all  problem 
categories  was  extremely  low  in  comparison  to  SBl. 

It  is  clear  that  the  simple  heuristics  embedded  in  Multi-PCP  are  too  weak  to  com¬ 
pete  with  the  extended  tabu  search  of  Taillard.  At  the  same  time,  Taillard’s  best 
solutions  were  obtained  by  running  the  tabu  search  procedure  3  to  4  times  for  a  large, 
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Figure  5:  Mean  %  deviation  from  best  solution  for  increasing  job/machine  ratio 
large  benchmark  problems 


Table  4:  Mean  and  standard  deviation  of  CPU  seconds  for  procedures,  Multi-PCP  and 
SBl,  performed  on  the  large  benchmark  problems 


Job  X 
Machine 

Multi-PCP 

SBl 

mean 

cr 

mean 

CT 

15  x  15 

1.20 

0.08 

5.10 

1.39 

20  X  15 

3.42 

0.33 

7.63 

1.06 

30  X  15 

11.90 

0.85 

14.68 

1.72 

50  X  15 

68.11 

7.12 

141.33 

104.84 

20  X  20 

3.73 

0.32 

15.64 

2.91 

30  X  20 

15.51 

0.77 

31.58 

4.03 

50  X  20 

94.90 

6.28 

165.83 

94.75 

100  X  20 

857.36 

38.43 

- 

- 

pre-determined  number  of  iterations  Given  the  additional  (but  incomplete)  computa¬ 
tional  details  provided  in  (Taillard,  1989)  on  incremental  solution  progress  and  account¬ 
ing  for  computing  hardware  differences,  it  would  appear  that  Multi-PCP  consistently 
obtains  its  solution  in  less  time  than  it  takes  Taillard’s  procedure  to  obtain  a  solution  of 
comparable  quality  (particularly  on  problems  with  low  job-to-machine  ratios).  Though 
this  computational  advantage  is  lowered  by  the  more  efficient  tabu  procedure  recently 
reported  in  (Nowicki  and  Smutnicki,  1994),  Multi-PCP  would  nonetheless  seem  to  offer 
leverage  as  a  solution  seeding  mechanism  to  such  iterative  search  procedures. 


5  More  Complicated  Problem  Formulations 


The  above  study  relates  the  performance  of  Multi-PCP  to  state-of-the-art  makespan 
minimization  procedures;  perhaps  somewhat  surprising,  it  shows  that  Multi-PCP’s  use 
of  a  CSP  scheduling  model  in  conjunction  with  fairly  simple  search  control  heuris¬ 
tics  leads  to  respectable  cost /effect! veness  (although  certainly  not  overwhelming  all 
previously  reported  results).  In  this  section,  we  consider  a  complementary  issue:  the 
broader  applicability  of  CSP-based  solution  approaches,  and  Multi-PCP  specifically, 
to  more  idiosyncratic  problem  formulations.  As  indicated  earlier,  real-world  applica¬ 
tions  are  often  complicated  by  additional  temporal  synchronization  and  resource  usage 
constraints,  and  solution  procedures  which  rely  on  problem  structure  that  is  pecu¬ 
liar  to  the  canonical  job  shop  problem  formulation  are  of  little  use  in  such  contexts. 
CSP  scheduling  models,  alternatively,  are  based  on  very  general  representational  as¬ 
sumptions  and  naturally  extend  to  accommodate  richer  problem  formulations.  We 

^Provisions  are  made  for  early  termination  upon  detection  of  an  optimal  solution.  As  indicated  in 
Appendix  A,  optimal  solutions  were  found  for  1  of  the  10  15x15  problems  and  for  9  of  the  10  problems 
in  each  of  the  15x50  and  20x100  categories.  In  all  other  problems,  there  was  no  early  termination  of 
the  search. 
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demonstrate  this  generality  below  by  applying  the  Multi-PCP  procedure  to  another, 
less-structured  makespan  minimization  problem;  the  multi-product  version  of  the  hoist 
scheduling  problem  studied  in  (Yih,  1994). 

5.1  The  Hoist  Scheduling  Problem 

The  hoist  scheduling  problem  finds  its  origin  in  printed  circuit  board  (PCB)  electro¬ 
plating  facilities.  In  brief,  a  set  J  of  jobs,  J  =  { Ji, ...,  Jn}  each  require  a  sequence  of 
chemical  baths,  which  take  place  within  a  set  M  of  m  chemical  tanks,  M  =  {1,  ...,m}. 
Execution  of  a  particular  chemical  bath  operation  Oi  requires  exclusive  use  of  tank 
rrii.  The  processing  time  of  any  Oi  required  for  a  job  j  is  not  rigidly  fixed;  instead 
there  is  a  designated  minimum  time,  p™*",  that  j  must  stay  in  the  tank  for  the  bath  to 
accomplish  its  intended  effect  and  a  maximum  time,  over  which  product  spoilage 

occurs.  All  jobs  move  through  the  chemical  tanks  in  the  same  order,  though  a  given 
job  may  require  only  a  subset  of  the  baths  and  thus  “skip”  processing  in  one  or  more 
tanks  along  the  way.  All  job  movement  through  the  facility  is  accomplished  via  a  sin¬ 
gle  material  handling  hoist,  which  is  capable  of  transporting  a  job  initially  into  the 
system  from  the  input  buffer,  from  tank  to  tank,  and  finally  out  of  the  system  into 
the  output  buffer.  H  can  grip  only  a  single  job  at  a  time,  moves  between  any  two 
adjacent  stations  (input  buffer,  tanks,  or  output  buffer)  at  constant  speed  5,  and  has 
constant  loading  and  unloading  speeds,  L  and  £/,  at  any  tank  or  buffer.  The  facility 
itself  has  no  internal  buffering  capability;  thus  jobs  must  be  moved  directly  from  one 
tank  to  the  next  once  they  have  entered  the  system.  The  objective  is  to  maximize 
facility  throughput  (or  equivalently  minimize  makespan)  subject  to  these  process  and 
resource  constraints. 

Most  previous  work  in  hoist  scheduling  has  considered  simplified  versions  of  this 
problem.  The  single-product,  hoist  scheduling  problem  has  received  the  most  attention. 
In  this  special  case,  the  problem  can  be  reduced  to  one  of  finding  a  minimum  length 
cycle  of  hoist  operations,  which  can  then  be  repeated  over  time;  several  algorithms 
for  generating  optimal  (or  near-optimal)  cyclic  schedules  have  been  reported  (Phillips 
and  Unger,  1976;  Lieberman  and  Turksen,  1981;  Shapiro  and  Nuttle,  1988;  Lei  and 
Wang,  1991;  Armstrong  et  ah,  1994).  In  (Yih  and  Thesen,  1991;  Yih  et  ah,  1993),  a 
hoist  scheduling  problem  involving  a  multi-product  facility  is  considered,  but  without 
permitting  variance  in  job  routings  (i.e.  no  tank  skipping).  To  our  best  knowledge,  only 
(Yih,  1994)  has  reported  procedures  for  solving  the  general  hoist  scheduling  problem 
defined  above. 
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5.2  Representation  as  a  Constraint  Graph 

The  GTCN  formalism  introduced  in  Section  2.1  requires  only  slight  extension  to  model 
the  hoist  scheduling  problem.  Let’s  first  consider  representation  of  the  process  con¬ 
straints  of  any  given  job  j.  For  each  individual  operation  Oi  in  j’s  process  sequence,  we 
define  three  constraint  graph  nodes:  two  time  points,  representing  Oi's  start  and  end 
points,  and  an  interval,  representing  Oi  itself.  A  given  pair  of  start  and  end  points  are 
related  to  its  corresponding  interval  through  use  of  the  qualitative  constraints  starts 
and  finishes  respectively.  The  values  ultimately  assigned  to  the  time  points  of  any  Oi 
in  the  constraint  graph  will  represent  Ohs  scheduled  start  and  finish  times. 

The  duration  of  a  given  operation  Oi  is  modeled  by  specifying  a  metric  constraint 
between  its  start  and  end  points.  There  are  two  cases,  corresponding  to  the  two  types 
of  operations  that  must  be  interleaved  to  process  any  given  job.  If  Oi  is  a  required  tank 
operation,  the  constraint  is  specified  to  enforce  minimum  and  maximum 

allowable  times  in  tank  rm.  If  alternatively,  Oi  is  a  hoist  (transport)  operation,  then 
the  constraint  is  instead  a  function  of  the  distance  to  be  traversed  and  hoist  loading, 
traveling  and  unloading  speeds.  The  constraint  is  specified  in  this  case, 

where  mt  =  L  +  s  *  [destinatioui  —  origiui)  +  U.  For  convenience,  we  assume  a 
correspondence  between  a  given  tank’s  index  and  its  location,  and  assign  indices  0  and 
m  -|-  1  to  the  system’s  input  and  output  buffers  respectively  to  complete  this  mapping. 

The  process  sequence  for  any  given  job  j  is  specified  by  temporally  relating  the 
intervals  defined  for  j’s  constituent  tank  and  transport  operations.  Since  the  end  of 
any  given  operation  Oi  must,  by  definition,  coincide  with  the  start  of  the  meets 

constraint  is  used  to  establish  this  linkage.  Oi  meets  Oi+i  implies  that  eti  =  Dis¬ 

junctive  constraints  on  resource  usage  are  specified  as  relation  sets  between  operations 
that  require  the  same  resource.  Again  there  are  two  cases.  For  tank  operations  Oi  and 
Oj  where  mi  =  m^,  the  constraint  {before,  after}  is  introduced.  Note  that  the  relation 
set  used  for  this  purpose  in  the  canonical  job  shop  problem,  {before- or-meets,  after-or- 
met-by},  is  not  correct  here,  since  it  is  physically  impossible  to  switch  between  tank 
operations  without  intermediate  loading  and  unloading.  Synchronization  of  competing 
hoist  operations  is  the  only  remaining  issue.  In  this  case,  however,  basic  qualitative 
relations  are  insufficient,  as  they  do  not  allow  us  to  account  for  the  setup  time  that 
may  be  required  to  position  the  hoist  at  the  loading  location. 

To  overcome  this  limitation,  we  extend  our  representation  of  qualitative  constraints 
to  optionally  include  a  metric  quantifier.  For  our  purposes  here,  it  is  sufficient  to  include 
only  the  following  two  extended  relations:  before-or-meets[lagtime]  and  after-or-met- 
by [lagtime],  where  lagtime  >  0  designates  a  minimum  metric  separation  between  the 
related  intervals.  Thus,  whereas  the  constraint  Oi  before-or-meets  Oj  implies  eti  ^  stj, 
the  extended  constraint  Oi  before-or-meets[hij]  Oj  implies  eti  +  hij  <  stj.  For  each  pair 
of  hoist  operations  Oi  and  Oj  belonging  to  different  jobs,  we  specify  the  constraint  Oi 
{before-or-meets[hij],  after-or-met-by[hjiJ}  Oj,  where  hij  =  s*  \destinationi  —  originj\ 
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Figure  6:  Example  1 
and  hji  =  s  *  \destinationj  —  origiriil. 

To  illustrate,  Figure  6  shows  a  simple  example  problem  (taken  from  (Yih,  1994)) 
involving  2  jobs  to  be  processed  in  a  3  tank  system.  Each  operation  is  displayed  below 
the  tank  that  is  required,  and  gives  its  minimum  and  maximum  processing  times.  The 
corresponding  constraint  graph  model  is  given  in  Figure  7  (with  qualitative  relations 
abbreviated  as  follows:  starts  (s),  finishes(f),  meets  (m),  before  (b),  after  (bi),  before- 
or-meets(bm),  after-or-met-by  (bmi)). 

5,3  PCP  Extensions 

Given  a  GTCN  model  of  a  hoist  problem,  we  can  distinguish  two  sets  of  disjunctive 
constraints  that  must  be  resolved  during  scheduling:  those  that  relate  competing  tank 
operations,  and  those  relating  pairs  of  hoist  operations.  In  both  cases,  the  set  of  possible 
values  that  might  be  selected/assigned  differs  from  the  {before- or-meets,  after-or-met- 
by}  relation  set  used  earlier  in  the  canonical  job  shop  model;  and  proper  treatment  of 
these  extended  relation  sets  necessitates  some  adjustment  to  the  base  PCP  constraint 
satisfaction  procedure  that  underlies  Multi-PCP. 

Recall  from  Section  2.2,  that  PCP  relies  on  shortest  path  lengths  as  a  basic  indicator 
of  sequencing  flexibility.  In  essence,  the  shortest  path  from  eti  to  stj  for  operations  Oi 
and  Oj,  designated  spij,  indicates  the  current  maximum  feasible  separation  between 
these  two  points.  In  the  case  of  a  {before-or-meets,  after-or-met-by}  relation  set,  this 
is  a  reliable  accounting  of  remaining  flexibility.  If  spij  >  0  then  Oi  before-or-meets 
Oj  is  still  feasible  and  vice  versa  (providing  the  basis  for  dominance  conditions);  if 
spij  <  spji  then  there  is  less  flexibility  remaining  to  sequence  Oi  before  Oj  than  there 
is  to  do  the  opposite.  However,  in  the  context  of  sequence-dependent  setups,  shortest 
path  lengths  provide  only  a  partial  and  distorted  estimate  of  sequencing  flexibility. 
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Figure  7:  The  constraint  graph  for  example  1 


To  solve  this  problem,  we  generalize  the  basic  measure  of  flexibility  in  PCP  to 
incorporate  sequence-dependent  lag  times.  Assume  hij  to  be  the  lag  time  required  if 
Oi  is  processed  before  Oj,  and  hji  be  the  lag  time  required  if  Oj  is  processed  before  Oi 
(i.e.,  the  constraint  specified  in  the  network  is  Oi  {^before-or-meetsfhij],  after-or-met- 
hyfhjij}  Oj).  We  revise  the  dominance  conditions  and  search  control  heuristics  specified 
in  Section  2.2  by  simply  substituting  the  extended  calculation  {spij  -  hij)  for  spij 
and,  likewise,  substituting  {spji  —  hji)  for  spji-  Note  that  these  revised  definitions  also 
accommodate  the  \^before,  after^  relation  set  required  for  synchronizing  tank  operations 
(in  which  case,  hij  and  hji  are  both  set  to  the  smallest  possible  temporal  increment), 
as  well  as  the  basic  {before-or-meets,  after-or-met-by}  relation  set  (where  hij,  hji  =  0). 


5.4  Performance  Results 

In  this  section,  we  adapt  the  Multi-PCP  procedure  defined  in  Section  3  to  incorporate 
the  extensions  discussed  above,  and  examine  its  performance  characteristics.  We  fol¬ 
low  the  same  general  experimental  design  of  (Yih,  1994),  and  similarly  consider  hoist 
scheduling  in  a  PCB  electroplating  facility  with  5  chemical  tanks.  All  problems  defined 
consisted  of  100  jobs,  each  with  randomly  generated  routings  and  tank  processing  time 
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constraints,  and  all  assumed  to  be  simultaneously  available.  Since  material  flow  is  uni¬ 
directional,  differences  in  job  routings  correspond  to  which  and  how  many  tanks  are 
skipped.  Experiments  were  conducted  to  evaluate  performance  along  two  dimensions 
relating  to  facility  constraints  and  operation:  first  as  function  of  the  relative  speed 
of  the  hoist  to  mean  tank  processing  time,  and  second  as  a  function  of  the  degree  of 
flexibility  provided  by  tank  processing  time  constraints.  To  calibrate  results,  prob¬ 
lems  were  also  solved  using  the  hoist  scheduling  procedure  previously  developed  by 
Yih  (Yih,  1994).  Both  procedures  were  implemented  in  C  and  run  on  a  Sun  SPARC 
10  workstation.  Before  describing  the  experiments  in  more  detail  and  presenting  re¬ 
sults,  we  briefly  summarize  Yih’s  approach  and  the  hoist  scheduling  implementation 
of  Multi-PCP. 


5.4.1  Yih’s  Approach 

In  (Yih,  1994),  two  heuristic  procedures  for  hoist  scheduling  are  proposed  and  evalu¬ 
ated.  First,  to  provide  a  baseline  for  comparison,  a  simple  “basic  algorithm”  is  speci¬ 
fied.  In  brief,  the  basic  algorithm  constructs  a  feasible  schedule  one  job  at  a  time,  by 
repeatedly  extending  a  partial  solution  to  include  the  schedule  of  an  additional  job. 
In  adding  any  given  job  into  the  schedule,  only  minimum  tank  processing  times  are 
considered.  The  job’s  first  tank  operation  is  tentatively  scheduled  to  start  as  early  as 
possible  and  subsequent  operations  in  the  routing  are  sequentially  checked  for  feasibil¬ 
ity  with  respect  to  resource  constraints.  If  a  constraint  conflict  is  detected  (due  either 
to  tank  or  hoist  unavailability),  the  first  operation  is  delayed  by  an  amount  sufficient  to 
avoid  the  conflict  and  the  feasibility  check  is  repeated.  The  process  terminates  as  soon 
as  a  conflict  free  start  time  for  the  job  is  found,  and  the  partial  solution  is  extended  to 
include  the  corresponding  schedule. 

The  second  proposed  procedure,  which  we  will  refer  to  below  as  the  “Yih94”  algo¬ 
rithm,  extends  the  basic  algorithm  to  exploit  any  scheduling  flexibility  that  is  provided 
by  maximum  tank  processing  time  constraints.  Specifically,  the  job  scheduling  cycle  is 
augmented  as  follows.  Upon  detection  of  a  hoist  availability  conflict  during  the  feasi¬ 
bility  checking  process,  a  backward  pass  is  first  made  through  the  routing  to  enlarge 
the  durations  of  previous  tank  operations  where  maximum  processing  time  constraints 
permit.  If  the  conflict  is  eliminated  by  this  action,  then  the  forward  feasibility  check 
is  resumed.  If  not,  the  job  start  time  is  delayed  as  before.  Experimental  results  re¬ 
ported  in  (Yih,  1994)  indicated  that  this  enhanced  algorithm  produced  schedules  with 
distinguishably  lower  makespan  than  those  generated  by  the  basic  algorithm  across  a 
variety  of  system  configurations.  Owing  to  its  simple,  myopic  nature,  the  Yih  algo¬ 
rithm  also  operates  very  efficiently;  100  job,  5  tank  problems  are  solved  in  seconds  with 
our  implementation. 
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5.4.2  A  Multi-PCP  Implementation  for  Hoist  Scheduling 

To  adapt  the  Multi-PCP  procedure  defined  in  Section  3  for  hoist  scheduling,  we  aug¬ 
ment  the  “Simple  PCP”  procedure  to  incorporate  the  extended  dominance  conditions 
and  search  control  heuristics  discussed  above,  and  specify  new  sub-procedures  for  deter¬ 
mining  lower  and  upper  bounds,  di  and  djj-^  on  the  common  deadline  interval  considered 
within  the  k-iteration  search.  To  establish  dx,,  we  simply  compute  the  minimum  total 
required  processing  time  (including  hoist  operations)  for  each  job  and  take  the  maxi¬ 
mum.  To  establish  du:  we  run  Yih’s  “basic  algorithm”  and  take  the  makespan  of  its 
result.  As  in  the  job  shop  scheduling  experiments,  the  number  of  iterations,  /c,  is  set 
to  8  and  not  varied. 

There  is  one  further  pragmatic  consideration.  Given  the  more  complex  structure  of 
the  constraint  graph  required  to  model  the  hoist  scheduling  problem,  it  is  not  possible 
to  validly  exploit  the  special-case,  “slack-based”  computations  of  (Erschler  et  ah,  1976) 
(see  Section  2.2)  as  a  surrogate  for  shortest-path  length  calculations  within  PCP  (as 
it  is,  for  example,  in  the  case  of  canonical  job  shop  problem  formulations).^  This 
significantly  raises  the  computational  overhead  of  constraint  propagation;  preliminary 
experiments  with  smaller,  50  job,  5  tank  problems  were  found  to  yield  run  times  of 
about  30  minutes. 

To  provide  a  more  computationally  competitive  alternative  to  Yih’s  algorithm,  a 
simple  problem  decomposition  method,  conceived  initially  by  (Ashour,  1967)  and  ex¬ 
tensively  studied  in  (Hirabayashi  et  al.,  1994),  is  introduced.  Under  this  scheme,  a 
solution  is  generated  by: 

1.  Partitioning  the  total  set  of  n  jobs  into  [n/b\  subproblems  with  >  6  jobs, 

2.  Solving  each  subproblem  independently  with  Multi-PCP  to  obtain  partial 

schedules,  and 

3.  Randomly  combining  the  [^/5J  partial  schedules  to  produce  an  n  job  schedule. 

In  the  experiments  reported  below,  a  subproblem  size  of  6  =  10  was  uniformly  adopted. 
For  these  100  job,  5  tank  problems,  overall  solution  time  was  about  100  seconds. 

^Interestingly,  the  need  for  explicit  shortest  path  length  computations  is  not  a  consequence  of  the 
extensions  introduced  to  accommodate  sequence-dependent  setups,  but  is  due  instead  the  presence 
of  bounded  interval  constraints  between  time  points,  i.e.,  minimum  and  maximum  processing  time 
constraints.  Slack-based  computations  remain  a  legitimate  computational  simplification,  for  example, 
in  solving  canonical  job  shop  problems  with  sequence-dependent  setups. 
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5.4.3  Sensitivity  to  Relative  Hoist  Speed 


One  dimension  along  which  PCB  Electroplating  system  configurations  can  vary  is  in 
the  relative  speeds  of  hoist  and  tank  operations.  Accordingly,  the  sensitivity  of  schedul¬ 
ing  performance  to  this  ratio  provides  one  interesting  basis  for  comparative  analysis. 
We  first  contrast  the  performance  of  Multi-PCP  and  Yih94  on  problem  sets  designed 
to  vary  the  ratio  7  =  where  p™*"  is  the  mean  minimum  processing  time  of 

tank  operations  and  s  is  the  speed  of  the  hoist  in  moving  between  adjacent  system 
locations.  A  total  of  100  problems  were  randomly  generated  according  to  the  following 
experimental  design: 

•  For  each  problem,  each  job  routing  was  generated  by  first  choosing  some  number 
of  tanks  to  skip  from  U[0,4],  and  then  selecting  (without  replacement)  which 
tanks  to  skip. 

•  10  different  problem  sets,  each  consisting  of  10  problems,  were  generated  by 
drawing  minimum  processing  times  for  tank  operations  from  the  following  10 
uniform  distributions:  U[5,15],  U[15,25],  U[25,35],  U[35,45],  U[45,55],  U[55,65], 
U[65,75],  U[75,85],  U[85,95],  and  U[95,105]. 

•  For  a  given  tank  operation  Oi,  was  defined  as  p  x  p™*”,  where  p  is  the 

tolerance  factor,  p  was  held  constant  at  1.0  in  this  experiment. 

•  For  all  problems,  we  assumed  a  hoist  travel  speed  s  of  2  and  hoist  loading  or 
unloading  time  of  1.5.  Thus,  the  10  problem  sets  reflect  7  ratios  of  5,  10,  15,  20, 
25,  30,  35,  40,  45,  and  50. 

To  characterize  results,  we  take  the  same  approach  as  in  (Yih,  1994);  we  use  the 
“basic  algorithm”  to  provide  a  benchmark  solution  and  compute  the  improvement  rate 
provided  by  both  Multi-PCP  and  Yih94.  The  precise  calculation  is: 

BasicM  ~  Xm 

I TTtpT' OZ) GTThGTli  icCL'tG’  —  ^ 

BasiCM 

where  Basicu  is  the  makespan  of  the  basic  algorithm’s  solution,  and  Xm  is  the 
makespan  of  the  solution  generated  by  procedure  Y,  AT  G  {Multi-PCP,  Yih94}. 

Figure  8  graphically  displays  the  average  results  obtained  for  each  problem  set  with 
Multi-PCP  and  Yih94,  ordered  by  increasing  7  ratio.  Both  procedures  are  seen  to  gen¬ 
erate  the  largest  improvement  for  values  of  7  in  the  range  of  [10,25],  with  improvement 
rates  degrading  as  7  becomes  larger  or  smaller.  In  the  case  of  Yih94,  no  improvement  is 
obtained  at  either  of  the  extreme  points  tested.  Multi-PCP,  alternatively,  yields  an  im¬ 
provement  rate  of  8%  at  the  smallest  7  value,  and  as  7  becomes  increasingly  larger,  its 
improvement  rate  stabilizes  at  about  15%.  Across  all  experiments,  Multi-PCP  is  seen 
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Figure  8:  Solution  improvement  rates  for  increasing  ratio  of  mean  processing  time  to 
hoist  speed 

to  produce  solutions  that,  on  average,  are  15%  better  (in  relation  to  the  benchmark 
solution)  than  those  obtained  with  Yih94. 

The  behavior  of  Yih94  at  extreme  7  ratios  is  predictable.  As  the  hoist  is  oper¬ 
ated  at  increasingly  slower  speeds  in  relation  to  tank  processing  time,  hoist  availability 
becomes  an  increasingly  dominant  constraint.  The  use  of  processing  time  flexibility 
becomes  increasingly  insufficient,  by  itself,  as  a  means  of  overcoming  hoist  delays,  and 
the  behavior  of  Yih94  increasingly  degenerates  to  that  of  the  basic  algorithm.  Alterna¬ 
tively,  at  increasingly  faster  relative  hoist  speeds,  hoist  availability  becomes  increasingly 
less  constraining.  Opportunities  to  exploit  processing  time  flexibility  correspondingly 
decrease,  again  leading  to  increased  reliance  on  the  basic  algorithm. 


5.4.4  Sensitivity  to  Duration  Flexibility 

To  consider  how  processing  time  flexibility  affects  scheduling  performance,  the  initial 
100-problem  experimental  design  was  repeated  for  three  additional  settings  of  p,  the 
tolerance  parameter  used  to  determine  for  any  generated  tank  operation  Oi-  Fig¬ 
ures  9  (a)  and  (b)  show  the  performance  results  obtained  with  Multi-PCP  and  Yih94, 
respectively,  for  problems  generated  with  p  settings  of  0.0  (no  flexibility),  0.5,  1.0  (the 
original  setting)  and  2.0. 

As  expected,  the  performance  of  both  procedures  is  seen  to  improve  as  the  duration 
flexibility  provided  by  minimum  and  maximum  processing  time  constraints  is  increased. 
Multi-PCP  does  a  better  job  of  exploiting  duration  flexibility  across  all  configurations. 
At  7  =  15,  the  difference  is  most  extreme.  With  this  system  configuration,  Multi-PCP 
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Figure  9:  Solution  improvement  rates  as  processing  time  flexibility  is  varied 
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obtained  a  30%  increase  in  solution  improvement  rate  as  p  is  raised  from  0  to  2  while  the 
corresponding  increase  in  improvement  rate  observed  for  Yih94  over  the  same  problem 
sets  was  just  6%.  Interestingly,  Multi-PCP  is  found  to  improve  the  benchmark  solution 
even  in  absence  of  processing  time  flexibility  (/>  =  0)  for  all  problem  sets  except  the 
extreme  7  =  5  problem  set,  where  hoist  availability  is  most  problematic;  in  this  case 
Multi-PCP’s  improvement  rate  increases  to  10%  as  hoist  availability  becomes  less  and 
less  constraining.  Yih94,  of  course,  yields  no  improvement  at  p  =  0,  since  it  degenerates 
to  the  basic  algorithm  when  there  is  no  processing  time  flexibility. 

More  generally,  the  results  obtained  in  these  experiments  illustrate  the  difficulty 
of  accommodating  complex,  interacting  constraints  within  scheduling  procedures  that 
proceed  via  explicit  commitment  to  specific  start  times.  In  the  case  of  the  Yih94 
algorithm,  a  priori  design  decisions  are  made  as  to  when  it  is  productive  to  extend 
processing  times  versus  delay  job  starts  and  what  specific  extend  or  delay  decisions 
should  be  taken.  Most  frequently,  however,  the  best  decision  actually  involves  some 
combination  of  these  two  alternatives.  This  tradeoff  is  not  considered  by  Yih94,  and, 
while  extensions  to  consider  additional  types  of  decisions  would  appear  possible,  it  is  not 
clear  how  one  could  extend  such  a  “fixed  times”  solution  approach  to  generally  address 
this  tradeoff  in  a  cost /effective  manner.  Multi-PCP,  alternatively,  does  not  suffer  from 
this  limitation.  By  operating  instead  in  the  space  of  sequencing  decisions,  this  tradeoff 
is  naturally  and  directly  considered.  There  is  no  need  to  design  the  algorithm  to  reason 
explicitly  about  the  types  of  constraints  involved;  decision-making  can  instead  be  based 
strictly  on  their  emergent  influence  on  the  evolving  partial  solution. 


6  Summary  and  Conclusion 


In  this  paper,  we  have  examined  the  applicability  of  constraint  satisfaction  problem 
solving  (CSP)  techniques  to  the  long  studied  problem  of  makespan  minimization,  argu¬ 
ing  their  utility  as  a  basis  for  the  design  of  cost /effective  scheduling  procedures  which 
naturally  accommodate  the  idiosyncratic  constraints  that  complicate  real-world  ap¬ 
plications.  We  reformulated  the  makespan  minimization  problem  as  a  search  for  the 
least  feasible  common  deadline  d,  and  proposed  a  simple  search  procedure  to  approx¬ 
imate  d  by  solving  a  series  of  deadline  scheduling  problems.  A  CSP  procedure  called 
PCP,  which  constructs  a  schedule  by  incrementally  adding  precedence  constraints  to  a 
temporal  constraint  graph  representing  the  input  problem,  was  embedded  as  the  core, 
deadline  scheduling  sub-procedure.  The  overall  procedure,  referred  to  as  Multi-PCP, 
was  then  empirically  evaluated  to  assess  its  performance  characteristics. 

To  calibrate  its  cost/effectiveness  as  an  approximate  solution  procedure,  the  perfor¬ 
mance  of  Multi-PCP  was  contrasted  with  that  of  recently  developed  shifting  bottleneck 
and  tabu  search  procedures  on  previously  studied  benchmark  problem  sets.  Despite  its 
fairly  simple  heuristics,  Multi-PCP  was  found  to  perform  quite  well.  On  the  Fisher  and 
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Thompson/ Lawrence  benchmarks,  Multi-PCP  produced  overall  results  comparable  to 
SB3  at  a  computational  cost  comparable  to  SBl  (Applegate’s  version).  Multi-PCP  was 
found  to  perform  best  on  problems  with  low  job-to-machine  ratios,  achieving  perfor¬ 
mance  comparable  to  SB4  on  the  10x10  problems.  On  the  larger  Taillard  benchmarks, 
Multi-PCP’s  sensitivity  to  job-to-machine  ratio  was  further  confirmed;  on  problems 
with  low  ratios,  Multi-PCP  consistently  produced  better  solutions  than  Applegate’s 
SBl  (the  only  shifting  bottleneck  procedure  we  were  able  to  compare),  while  the  in¬ 
verse  was  observed  at  higher  ratios  (for  those  problems  that  the  SBl  implementation 
was  able  to  successfully  solve).  In  no  cases,  did  Multi-PCP  (or  SBl)  produce  so¬ 
lutions  as  good  as  those  obtained  with  Taillard’s  tabu  search  procedure.  However, 
Multi-PCP  solutions  were  obtained  in  less  time  than  it  took  to  generate  comparable 
solutions  within  the  tabu  search,  suggesting  its  potential  as  a  seeding  mechanism  for 
such  extended  search  procedures. 

To  demonstrate  the  generality  of  the  approach,  we  also  considered  application  of 
Multi-PCP  to  the  multi-product,  hoist  scheduling  problem  previously  studied  by  Yih. 
This  problem  requires  satisfaction  of  a  more  complicated  set  of  constraints  than  does 
the  canonical  job  shop  scheduling  formulation  of  the  benchmark  problems,  including 
simultaneity  in  the  starts  and  ends  of  interleaved  material  transport  and  manufacturing 
operations,  non-rigid  processing  times,  and  sequence- dependent  hoist  travel  times.  It  is 
in  this  problem  context  that  the  power  of  the  underlying  constraint  satisfaction  model 
is  most  evident.  With  only  minor  extension  to  account  for  sequence-dependent  setups, 
Multi-PCP  was  shown  to  be  directly  applicable  to  this  more  idiosyncratic  makespan 
minimization  problem.  Though  computational  considerations  necessitated  use  of  a 
problem  decomposition  scheme  to  achieve  realistic  problem  scale,  experimental  re¬ 
sults  nonetheless  indicated  uniform,  sizable  improvement  in  solution  quality  over  other 
known  solution  procedures  across  a  wide  range  of  system  configurations. 
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Table  5:  Makespan  results  for  small  benchmark  problem  set 


Problem 

Job  X 
Machine 

Optimal 

Value 

Multi- 

PCP 

Multi-PCP 
w/  Relax 

SBl 

SB3 

SB4 

mt06 

6x6 

55 

55 

55 

59 

55 

55 

mtlO 

10  X  10 

930 

949 

949 

952 

981 

940 

mt20 

20  X  5 

1165 

1267 

1192 

1228 

1199 

1199 

lal 

666 

666 

666 

666 

666 

666 

la2 

655 

670 

670 

684 

667 

667 

la3 

10  X  5 

597 

619 

617 

605 

626 

626 

la4 

590 

600 

600 

603 

593 

593 

la5 

593 

593 

593 

593 

593 

593 

la6 

926 

926 

926 

926 

926 

926 

la7 

890 

890 

890 

890 

890 

890 

laS 

15  X  5 

863 

863 

863 

863 

863 

863 

la9 

951 

951 

951 

951 

951 

951 

lalO 

958 

958 

958 

958 

958 

958 

lall 

1222 

1222 

1222 

1222 

1222 

1222 

lal2 

1039 

1039 

1039 

1039 

1039 

1039 

lal3 

20  X  5 

1150 

1150 

1150 

1150 

1150 

1150 

lal4 

1292 

1292 

1292 

1292 

1292 

1292 

lal5 

1207 

1207 

1207 

1207 

1207 

1207 

lal6 

945 

982 

982 

1076 

961 

961 

lal7 

784 

787 

787 

829 

796 

796 

lal8 

10  X  10 

848 

886 

886 

855 

866 

861 

lal9 

842 

852 

852 

863 

902 

878 

la20 

902 

922 

922 

918 

932 

922 

927 

984 

965 

968 

954 

954 

la23 

15  X  10 

1032 

1041 

1041 

1071 

1032 

1032 

la24 

935 

983 

983 

1015 

976 

976 

la25 

977 

1026 

1026 

1061 

1012 

1012 

la26 

1218 

1272 

1272 

1393 

1239 

1224 

la28 

1216 

1275 

1275 

1281 

1245 

1245 

la30 

1355 

1356 

1356 

1355 

1355 

1355 

la31 

1784 

1784 

1784 

1784 

1784 

1784 

la32 

1850 

1850 

1850 

1850 

1850 

1850 

la33 

1719 

1722 

1722 

1719 

1719 

1719 

la34 

1721 

1744 

1744 

1721 

1721 

1721 

la35 

1888 

1888 

1888 

1888 

1888 

1888 

la36 

1268 

1321 

1321 

1326 

1319 

1319 

la37 

1397 

1446 

1466 

1471 

1425 

1425 

la39 

mm 

1233 

1286 

1286 

1301 

1278 

1278 

■■ 

1233 

1272 

1272 

1347 

1266 

1262 
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Table  6:  %  deviation  from  optimal  solution  for  small  benchmark  problem  set 


Problem 

Job  X 
Machine 

Multi- 

PCP 

Multi-PCP 
w/  Relax 

SBl 

SB3 

SB4 

mt06 

6x6 

0 

0 

7.27 

0 

0 

mtlO 

10  X  10 

2.04 

2.04 

2.37 

5.48 

1.08 

mt20 

20  X  5 

8.76 

2.32 

5.41 

2.92 

2.92 

lal 

0 

0 

0 

0 

0 

la2 

2.29 

2.29 

4.43 

1.83 

1.83 

la3 

10  X  5 

3.69 

3.35 

1.34 

4.86 

4.86 

la4 

1.69 

1.69 

2.20 

0.51 

0.51 

la5 

0 

0 

0 

0 

0 

la6 

0 

0 

0 

0 

0 

la7 

0 

0 

0 

0 

0 

la8 

15  X  5 

0 

0 

0 

0 

0 

la9 

0 

0 

0 

0 

0 

lalO 

0 

0 

0 

0 

0 

lall 

0 

0 

0 

0 

0 

lal2 

0 

0 

0 

0 

0 

lal3 

20  X  5 

0 

0 

0 

0 

0 

lal4 

0 

0 

0 

0 

0 

lal5 

0 

0 

0 

0 

0 

lal6 

3.92 

3.92 

13.86 

1.69 

1.69 

lal7 

0.38 

0.38 

5.74 

1.53 

1.53 

lal8 

10  X  10 

4.48 

4.48 

0.83 

2.12 

1.53 

lal9 

1.19 

1.19 

2.49 

7.13 

4.28 

la20 

2.22 

2.22 

1.77 

3.33 

2.22 

la22 

6.15 

4.10 

4.42 

2.91 

2.91 

la23 

15  X  10 

0.87 

0.87 

3.78 

0 

0 

la24 

5.13 

5.13 

8.56 

4.39 

4.39 

la25 

5.02 

5.02 

8.60 

3.58 

3.58 

la26 

4.43 

4.43 

14.37 

1.72 

0.49 

la28 

20  X  10 

4.85 

4.85 

5.35 

2.38 

2.38 

la30 

0.07 

0.07 

0 

0 

0 

la31 

0 

0 

0 

0 

0 

la32 

0 

0 

0 

0 

0 

la33 

30  X  10 

0.17 

0.17 

0 

0 

0 

la34 

1.34 

1.34 

0 

0 

0 

la35 

0 

0 

0 

0 

0 

la36 

4.18 

4.18 

4.57 

4.02 

4.02 

la37 

15  X  15 

4.94 

4.49 

5.30 

2.00 

2.00 

la39 

4.30 

4.30 

5.52 

3.65 

3.65 

la40 

3.16 

3.16 

9.25 

2.68 

2.35 

37 


Table  7:  Computation  times  (in  seconds)  for  small  benchmark  problem  set 


Problem 

Job  X 
Machine 

Multi- 

PCP 

SBl 

SB3 

SB4 

mt06 

6x6 

1.45 

mtlO 

10  X  10 

0.38 

2.21 

7.76 

mt20 

20  X  5 

1.38 

1.99 

3.62 

lal 

0.05 

0.12 

0.13 

la2 

0.23 

0.13 

0.51 

0.81 

la3 

10  X  5 

0.17 

0.10 

0.43 

0.68 

la4 

0.19 

0.11 

0.37 

0.63 

la5 

0.02 

0.13 

0.61 

0.61 

la6 

0.03 

0.15 

la7 

0.53 

0.17 

0.11 

0.11 

laS 

15  X  5 

0.48 

0.18 

0.20 

la9 

0.05 

0.12 

0.11 

0.11 

lalO 

0.02 

0.15 

0.10 

lall 

0.03 

0.18 

0.19 

0.19 

lal2 

0.03 

0.12 

0.12 

0.12 

lal3 

20  X  5 

0.05 

0.18 

0.13 

0.13 

lal4 

0.03 

0.13 

0.09 

lal5 

0.20 

0.20 

0.18 

0.18 

lal6 

0.27 

0.48 

1.34 

3.59 

lal7 

0.25 

0.65 

1.57 

2.97 

lal8 

10  X  10 

0.23 

0.62 

1.74 

lal9 

0.30 

0.78 

1.67 

la20 

0.23 

0.80 

1.71 

la22 

1.15 

3.59 

6.52 

la23 

15  X  10 

1.10 

1.28 

1.99 

1.99 

la24 

0.95 

DEI 

3.92 

7.19 

la25 

0.97 

tm 

4.15 

7.25 

la26 

2.55 

1.33 

4.55 

10.77 

la28 

2.37 

1.83 

5.71 

10.41 

la30 

2.62 

1.67 

1.33 

1.33 

la31 

5.77 

3.04 

3.45 

3.45 

la32 

0.12 

2.38 

2.74 

2.74 

la33 

30  X  10 

8.42 

1.67 

4.63 

4.63 

la34 

8.63 

3.04 

6.99 

6.99 

la35 

1.30 

2,77 

2.86 

2.86 

la36 

1.08 

4.07 

14.01 

28.14 

la37 

15  X  15 

1.23 

2.93 

13.22 

26.84 

la39 

1.47 

5.78 

12.84 

25.59 

la40 

1.37 

5.40 

13.30 

26.58 

Table  8:  Makespan  results  for  large  benchmark  problem  set 


Job  X 

Best 

Multi- 

Multi-PCP 

Job  X 

1  Best 

Multi- 

Multi-PCP 

Machine 

Value 

PCP 

w/  Relax 

SBl 

Machine 

:  Value 

PCP 

w/  Relax 

1231* 

1280 

1280 

1360 

2064 

2323 

2323 

1252 

1322 

1308 

1367 

1983 

2205 

2205 

1223 

1288 

1288 

1289 

1896 

2165 

2165 

1181 

1260 

1253 

1289 

2031 

2269 

2252 

15  X  15 

1234 

1306 

1306 

1359 

30  X  20 

2032 

2215 

2215 

1243 

1320 

1320 

1314 

2057 

2269 

2269 

1228 

1307 

1307 

1322 

1947 

2104 

2104 

1221 

1289 

1289 

1345 

2005 

2277 

2277 

1289 

1366 

1366 

1437 

2013 

2337 

2330 

1261 

1334 

1334 

1331 

1973 

2188 

2188 

1376 

1522 

1488 

1481 

2760* 

2948 

2948 

1381 

1495 

1495 

1503 

2756* 

2976 

2964 

1367 

1478 

1478 

1521 

2717* 

2898 

2891 

1355 

1452 

1452 

1540 

2839* 

2905 

2905 

20  X  15 

1366 

1486 

1486 

1532 

50  X  15 

2689 

2944 

2944 

1371 

1478 

1478 

1511 

2781* 

3023 

3023 

1480 

1609 

1609 

1605 

2943* 

3178 

3156 

1432 

1503 

1503 

1532 

2885* 

3131 

3131 

1361 

1440 

1440 

1504 

2655* 

2943 

2938 

1373 

1441 

1441 

1535 

i 

2723* 

2904 

2904 

1663 

1764 

1764 

1814 

2921 

3177 

3177 

1626 

1739 

1739 

1776 

3002 

3311 

3244 

1574 

1729 

1687 

1822 

2835 

3059 

3047 

1660 

1773 

1773 

1822 

2775 

2982 

2951 

20x  20 

1598 

1742 

1742 

1847 

50  X  20 

2800 

3061 

3061 

1679 

1759 

1759 

1861 

2914 

3188 

3187 

1704 

1862 

1862 

1857 

2895 

3165 

3165 

1626 

1749 

1749 

1809 

2835 

3048 

3048 

1635 

1761 

1761 

1771 

3097 

3419 

3419 

1614 

1743 

1743 

1729 

3075 

3349 

3349 

1770 

2015 

2001 

2017 

5464* 

5804 

5804 

1853 

2033 

2019 

2061 

5181* 

5378 

5378 

1855 

2053 

2053 

1987 

5568* 

5885 

5885 

1851 

2071 

2071 

1996 

5339* 

5553 

5553 

30  X  15 

2007 

2106 

2106 

1965 

100  X  20 

5392* 

5743 

5743 

1844 

2031 

2016 

1965 

5342* 

5733 

5733 

1822 

1954 

1954 

1976 

5436* 

5621 

5621 

1714 

1853 

1853 

1915 

5394* 

5618 

5618 

1824 

2010 

2010 

1890 

5358* 

5563 

5563 

1723 

1931 

1931 

1838 

5213 

5413 

5413 

These  indicate  optimal  values. 


Table  9:  %  deviation  from  best  solution  for  large  benchmark  problem  set 


Job  X 
Machine 

Multi- 

PCP 

Multi-PCP 
w/  Relax 

SBl 

Job  X  i 
Machine  I 

Multi- 

PCP 

Multi-PCP 
w/  Relax 

SBl 

3.98 

3.98 

10.48 

12.55 

12.55 

13.52 

5.59 

4.47 

9.19 

11.20 

11.20 

10.89 

5.31 

5.31 

10.55 

14.19 

14.19 

11.97 

6.69 

6.10 

9.14 

11.72 

10.88 

17.77 

15  X  15 

5.83 

5.83 

10.13 

30x20 

9.01 

9.01 

11.32 

Table  10:  Computation  times  (in  seconds)  for  large  benchmark  problem  set 


Job  X 

Multi- 

Job  X 

Multi- 

Machine 

PCP 

SBl 

Machine 

PCP 

SBl 

1.15 

4.08 

14.65 

28.32 

1.18 

5.88 

14.05 

37.95 

1.17 

6.48 

15.67 

35.42 

1.13 

4.52 

16.32 

28.97 

15  X  15 

1.33 

3.22 

30  X  20 

14.85 

34.05 

1.20 

5.05 

16.52 

32.73 

1.17 

3.37 

15.87 

34.70 

1.27 

6.88 

15.60 

30.47 

1.33 

6.93 

15.68 

28.30 

1.10 

4.55 

15.90 

24.92 

3.53 

7.98 

73.88 

269.91 

3.05 

9.62 

72.20 

139-21 

3.15 

6.65 

66.68 

316.89 

3.18 

6.20 

72.46 

47.40 

20  X  15 

3.67 

8.42 

50  X  15 

64.68 

- 

3.32 

6.97 

50.43 

46.26 

3.63 

7.77 

70.48 

37.30 

4.10 

7.77 

75.31 

175.06 

3.54 

8.43 

68.66 

49.80 

3.06 

6.48 

66.28 

190.11 

3.40 

11.62 

93.43 

125.86 

4.32 

18.48 

91.43 

150.18 

3.40 

15.92 

89.01 

251.82 

3.77 

13.90 

90.74 

98.21 

20  X  20 

3.78 

13.75 

50  X  20 

103.90 

382.80 

3.80 

13.78 

105.43 

90.95 

4.05 

13.52 

95.35 

97.71 

3.77 

15.53 

100.61 

145.39 

3.22 

19.73 

91.33 

- 

3.83 

20.20 

87.72 

149.58 

12.28 

12.45 

944.45 

“ 

12.54 

15.97 

832.63 

- 

10.98 

14.32 

870.70 

- 

11.95 

12.11 

798.67 

- 

30  X  15 

10.00 

14.17 

100  X  20 

865.13 

- 

12.59 

16.20 

846.10 

- 

11.60 

15.72 

836.57 

- 

12.77 

17.68 

864.87 

- 

12.43 

14.15 

877.18 

- 

11.82 

14.02 

837.33 

- 

41 


